Testing for topological order in variational wavefunctions for Z2 spin liquids 
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We determine the conditions under which a spin- liquid Mott insulator |0) defined by a Gutzwiller 
projected BCS state at half-filling is fractionalized. We construct a trial vison (Z^ vortex) 
state \V) by projecting an hc/2e vortex threading the hole of a cylinder/torus and examine its 
overlap with |0) using analytical and numerical calculations. We find that generically the overlap 
vanishes in the thermodynamic limit, so the spin-liquid is Z2 fractionalized. We point out the 
relevance of these results to numerical studies of Hubbard-like models and spin models which have 
been recently reported to possess spin liquid phases. We also consider possible implications for 
flux-trapping experiments that have tested for Z2 fractionalization in underdoped high temperature 
superconductors. 
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I. INTRODUCTION 

Spin liquid Mott insulators have long been viewed as 
candidates for the ground state of certain frustrated mag- 
netic systems^ in dimensions D > 2. They have been of 
interest in the context of high temperature superconduc- 
tivity following Anderson's proposal that doping such 
insulators may be a novel route to superconductivity. 
Progress in our theoretical understanding of a certain 
class of spin liquids has come from two lines of attack 
over the past decade. First, working from the point of 
effective theories, it has been shown that spin-liquid in- 
sulators in dimensions D > 2 may emerge as deconfincd 
quantum phases of Z2 gauge theories 3 ^^ at zero temper- 
ature. In this phase, an S = 1 excitation, which is the el- 
ementary excitation of a conventional magnet, may break 
up into two 5 = 1/2 particles, called spinous, which are 
minimally coupled to an Ising gauge field; hence also the 
term u Zz fractionalized spin liquids" for such insulators. 
The more recent line of progress has been in the construc- 
tion of simple model Hamiltonians which have no special 
gauge symmetries but which can nevertheless be shown 
to possess such Z2 fractionalized phases — these include 
quantum dimer models on the triangular 6 and Kagome 7 
lattices, and certain Hubbard-like boson models^. 

There are a few experimental indications of spin liquid 
states in two dimensions. Recent low temperature NMR 
and susceptibility measurements^ in a quasi-2D organic 
insulator indicate no magnetic phase transitions down to 
AOmK. Similarly, no ordering transitions are found in 
Helium-3 adsorbed on graphite^ even at very low tem- 
peratures around 10/iK. These temperatures are nearly 
two orders of magnitude smaller than the estimated mag- 
netic exchange energy scale in these respective systems. 
Further, experimental signatures of spinons have been re- 
ported in neutron scattering experiments^ on CS2CUCI4. 

In the past few years, there has also been a lot of nu- 
merical work on microscopic models, which are not an- 



alytically tractable, searching for spin-liquid phases. In 
particular, exact diagonalization studies of a multiple- 
spin exchange model on a triangular latticed, and Monte 
Carlo studies of certain two-dimensional Hubbard-like 
modelsi^i^ indicate that such models may possess, in a 
regime of parameters, insulating phases with no obvious 
broken translational or broken spin rotational symme- 
tries. 

What kind of spin liquids could these experiments and 
numerical studies be probing? In order to show how to 
answer this question in the context of numerical studies, 
we focus on "topological order" which is a sharp test of 
whether the phases being accessed in the numerics are 
Z 2 fractionalized spin liquids. Topological order refers to 
degeneracies in the spectrum of the Hamiltonian which 
depend on the topology of the manifold on which the 
system lives. This may be understood as follows: if the 
system is in the deconfined phase of a Z2 gauge theory, it 
necessarily implies a novel gapped excitation in the bulk 
of the systen4*^*i&, namely the Ising vortex of the Z2 
gauge field called the "vison" . However, a vison thread- 
ing noncontractible loops of the manifold, i.e. threading 
the holes of the cylinder or torus, does not cost any en- 
ergy in the thermodynamic limit. This leads to a topolog- 
ical degeneracy of the low energy eigenstate a 4 i 5 ' 16 , which 
may be viewed as states with/ without visons threading 
the holes of the cylinder or torus. In summary, for a Z2 
fractionalized spin-liquid insulator we expect a two-fold 
degeneracy of the low energy spectrum on a cylinder and 
a four- fold degeneracy on a torus. If indeed the above nu- 
merical studies are probing Z2 fractionalized spin-liquid 
phases, it would be useful to check if they possess topo- 
logical order consistent with this phased 

In this paper, we address the following issues: (a) How 
can one think about the vison in terms of electronic co- 
ordinates? (b) Under what conditions is the vison well- 
defined, leading to ^-fractionalization and topological 
degeneracy in spin-liquid insulators? We show how one 
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may obtain an estimate of the length scale beyond which 
fractionalization is apparent, consider its implications for 
the cuprate superconductors, and also point out how one 
can use the method described here to test for topo- 
logical order in other numerical studies of Hubbard-like 
models, superconductors. Our work builds upon and ex- 
tends some of the results of Ivanov and Senthil 18 . 

Throughout this paper, we will discuss the above is- 
sues in the context of a specific spin-liquid state, namely 
a Gutzwiller projected d-wave BCS wavefunction on a 
square lattice. This wavefunction has been shown to 
provide a remarkably good description of the supercon- 
ducting state of the high Tc cuprates over a wide range 
of dopingiS, and the present study was carried out to 
examine its implications for the undoped insulator and 
explore connections with phenomenological theories of 
the cuprates involving Z2 fractionalization. However, our 
study serves more generally to illustrate the method of 
identifying Z 2 spin-liquid states either using Gutzwiller 
projected variational wavefunctions, or even other nu- 
merical methods as we discuss in Section 1X1 

II. CONSTRUCTING THE VISON 
WAVEFUNCTION 

We begin with a discussion of how to construct a vison 
wavefunction in terms of microscopic electronic degrees of 
freedom, which form the basis of the numerical studies. 
For conceptual clarity, consider first a square lattice of 
L 2 sites "wrapped" into a cylinder along the x direction. 
The spin- liquid ground states |0) of interest to us are 
given by the Gutzwiller projection of ./V-particle d-wave 
BCS states at half-filling (N = L 2 ): 

|0) = V\BCS) = 7>[^>(r - r')4 T 4,J W/2 |vac>. (1) 

r,r' 

Here the pair wavefunction is given explicitly by 

<p{p) = L- 2 exp( l k • p) [A k /(a + + A*)] (2) 

k 

with k's are chosen consistent with periodic boundary 
condition (PBC) along x. The pair wavefunction is 
parametrized in terms of two variational parameters- 8 
/1 and A which determine £k = e(k) — /1, with 
e(k) = — 2t (cos k x + cos k y ) — At' cos k x cos k y and Ak 
A (cosfej. — cos k y ) /2. The Gutzwiller projection opera- 
tor V = rir(l ~ n rt n ri) restricts the configurations to 
have exactly one electron per site. 

While the state before projection describes a super- 
conductor, the projected state is an insulator with one 
electron (or equivalently one spin) per site, and no long 
range magnetic order—. 

To test for Z2 fractionalization, we would like to con- 
struct a trial state for a vison threading a hole of the 
cylinder to check for topological order. We do this in two 
steps. First, we construct a wavefunction for an hc/2e 



vortex threading the hole of the cylinder. Next, we we 
Gutzwiller project this hc/2e vortex state and present 
four arguments for why gives us a good candidate for a 
vison state. 

To construct a superconducting state with one hc/2e 
vortex threading the hole of the cylinder, we modify the 
pairing function in the BCS ground state as 

<p(r - r') -> exp[iQ • (r + r')/2]^(r - r'), (3) 

where Q = 2-kx/L. This gives a center of mass mo- 
mentum Q to each Cooper pair and twists the phase of 
the condensate by 2tt going around the cylinder. How- 
ever, the requirement of a single-valued wavefunction 
when even one of the electrons of the Cooper pair (r 
or r') goes around the cylinder, constrains k's to sat- 
isfy aniiperiodic boundary conditions (APBC) along x 
in the Fourier transform in defining (Pa(t — r') via equa- 
tions analogous to Eas. l|ll2|l : hence the subscript A on 
<Pa(t — r ') m Eq. [3] above. Clearly the state with one 
hc/2e vortex has a current flow and it can also be writ- 
ten as 

\hc/2e) y = exp (i — ^— [n T (r) + n L (r)}) 

r 

xtVJ^^-rOc^c^J^lvac). (4) 

r,r' 

Gutzwiller projecting this state at half-filling fixes 
n|(r) + nj,(r) = 1, and the prefactor multiplying the 
wavefunction in Eq. Q then drops out as a trivial phase 
factor since it is independent of the spin configuration. 
This gives 

\V V ) = V\hc/2e) y =V[J2^P A (r- r')4 T 4'J N/ * (5) 

r,r' 

Thus the Gutzwiller projected hc/2e vortex state \V y ) 
simplifies to a form identical to Eq. Q with <p — > tpA, 
and corresponds to imposing antiperiodic boundary con- 
ditions on the electrons before projecting them. 

Why is this a good candidate for a vison state? First, a 
projected hc/e vortex with Q = Anx/L does not require 
any change in BC's on tp to ensure single- valuedness, and 
the phase factors again drop out at half-filling. It is thus 
identical to the ground state |0). Hence we only have 
two kinds of vortex states upon projection: all super- 
conducting states with an even number of hc/2e vortices 
threading the cylinder collapse onto the ground state |0), 
those with an odd number of hc/2e vortices threading 
the cylinder collapse onto \V y ). This renders manifest 
the Zi character of a projected hc/2e vortex — two vor- 
tices is the same as none — which makes them attractive 
candidates for visons. 

Second, we can look at the crystal momentum of the 
projected vortex. The vortex state before projection has 
a crystal momentum Q per singlet pair, or a total crys- 
tal momentum of QA^/2. Since the Gutzwiller projec- 
tion operator commutes with the translation operator, 
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the crystal momentum is unchanged by projection. Thus 
for a system with dimensions L x x L y , the projected vor- 
tex has a total crystal momentum P x — TrL y . Since the 
crystal momentum is only defined modulo 2tt, this means 
the projected hc/2e vortex wavefunction carries momen- 
tum 7r for odd- L y and no momentum for even-Ly. As 
shown elsewhere 21 *, this is precisely the momentum of a 
vison in an insulator described by a Z2 gauge theory with 
one spin per site. Thus the momentum quantum number 
of the projected hc/2e vortex is consistent with that of 
the vison. As discussed below, we will always work with 
even-L y in this paper. 

Third, one can consider a specific limit of the pair 
wavefunction, namely /1 — > —00, which corresponds to 
having singlet pairs only on neighboring sites, like in 
a nearest neighbor dimer model. The state |0) corre- 
sponds to a superposition of nearest neighbor dimer (sin- 
glet) configurations. In this limit, due to the antiperiodic 
boundary conditions on (^(r — r'), it is straightforward 
to show that every configuration in which an odd number 
of singlets lie on links such that one spin of the singlet is 
at a site (L x ,y) and the other spin lies at (1,2/) (taking 
into account all allowed y) has an amplitude which dif- 
fers in sign from its amplitude in the state |0). This is 
however identical to the well-known construction of the 
Z2 vison threading a cylinder in the context of dimer 
models 1 ^. 

Finally, in terms of a dual theory of vortices, the way 
to obtain a Z2 fractionalized insulator starting from a 
superconducting state is by condensing pairs of hc/2e 
vortices 2 ^, thus quantum disordering the superconduc- 
tor. In this case, the vison is the remnant of the hc/2e 
vortex which has not condensed. Since Gutzwillcr pro- 
jecting a superconductor at half-filling fixes the number 
of particles, it amounts to phase disordering the super- 
conductor, and this also leads one to suspect that the 
Gutzwiller projected hc/2e vortex is the vison. 

Generalizing the above discussion to a torus, |0) cor- 
responds to PBC along both x and y, \V X ) and | V^) cor- 
respond to PBC/APBC and \V xy ) to APBC along both 
directions. These proposed vison wavefunctions are de- 
signed to investigate the existence of disconnected topo- 
logical sectors as evidence for fractionalization as we dis- 
cuss in the following section. A vison excitation in the 
bulk, on the other hand, would need to have a Z2 flux 
piercing the plane of the lattice and would require careful 
consideration of the superconducting vortex core prior to 
projection, which we avoid here. 



III. HOW DO WE TEST FOR TOPOLOGICAL 
ORDER? 

A SC state with a vortex threading the hole of the 
cylinder/torus carries current and is orthogonal to the 
ground state. However, Gutzwiller projecting the BCS 
wavefunction at half-filling destroys SC order— and re- 
sults in an translationally invariant insulator—. The key 



question then is: does any remnance of vorticity sur- 
vive projection? If it does, then (V a \0) = (V a \Vp) = , 
for afi—x, y, xy, the vison is well-defined leading to topo- 
logical degeneracy^, and one is in a Z2 fractionalized 
phase. If, on the other hand, the proposed vison state 
has nonzero overlap with the state |0) in the thermo- 
dynamic limit, one is in a conventional unfractionalized 
phase. 

Since we have shown above that the vison state on 
a cylinder with odd-Ly carries crystal momentum 7r, it 
is trivially orthogonal to the ground state |0) even on 
a finite system. However this by itself does not tell us 
anything about topological order since if this state be- 
comes degenerate with the ground state it may mix with 
the ground state in the thermodynamic limit and break 
translational symmetry — for instance, this is what is 
known to happen in any gapped phase of an odd number 
of coupled spin chains (an "odd-leg ladder"), as is well- 
known from an extension 22 of the Lieb-Schultz-Mattis 
theorem in one dimension— However, in the limit of in- 
finite number of coupled chains, the order parameter in 
the broken symmetry state may survive leading to a con- 
ventional broken symmetry state in two dimensions, or 
it may vanish leading to a translationally invariant two- 
dimensional spin liquid 24 . These two possibilities cannot 
be distinguished by an overlap calculation with odd-L y . 

For purposes of testing for topological order, we are 
therefore really interested in the case of even-Ly where 
the trial vison state carries no crystal momentum, and it 
is thus not trivially orthogonal to |0) on finite-size sys- 
tems. Showing that the overlap between the ground state 
and the proposed vison state vanishes in the thermody- 
namic limit on cylinders or torii with even length in each 
direction is therefore a crucial step in establishing topo- 
logical order. 



IV. "PHASE DIAGRAM" AND SYMMETRIES 

To determine the conditions under which the vison sur- 
vives, it is useful to consider the full parameter space for 
projected BCS states at half-filling, which is the (i',/i) 
plane with t = 1 and A held fixed at some nonzero value 
in order to describe a RVB liquid of singlet pairs. Study- 
ing this general class of states will be useful since one can 
check if there are states which are not Z2 fractionalized 
in some regime of parameters, and ask how the transi- 
tion from a fractionalized regime to an unfractionalized 
regime is reflected in the overlap of the wavefunctions 
(V a \V p ). 

We use symmetry arguments to show that not all 
states in this space are distinct upon projection and 
|0(t',/x)) = |0(—t', —//)). In brief, we can change t' — > 
— fi — > —fi by a global particle-hole transformation 
c va — > (— l) x+y c) ra in the wavefunction, also redefining 
the I vac) since empty sites transform into doubly occu- 
pied ones. However, with exactly one particle per site 
such a particle-hole transformation interchanges T-spins 
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FIG. 1: (A) "Phase diagram" showing that |0) is Z2 frac- 
tionalized in the entire (t' , /i)-plane except at the bipartite 
symmetric point (0,0) (see text for details). (1) Shaded 
and unshaded half-planes are related by symmetry (t',/f) — > 
(— t', —fJ>)', (2) hatched region for large negative /i is the short- 
range RVB (resonating valence bond) limit; (3) analytical ar- 
guments for fi = ^bcs (schematic line) are described in the 
text; (4) crosses indicate some points where we numerically 
compute overlaps; (5) shaded ellipse around (—0.25,0) is the 
region relevant for high Tc SCs. (B) Energy of |0) for the Hub- 
bard model at half-filling for parameters relevant to cuprates 
(U/t=12,t' /t=-0.25 and optimal A=1.25t) versus fi/t. The 
optimal value is —0.3 < fj,/t ^ 0.3 within statistical errors. 

and j,-spins at every site for any configuration. Since 
the BCS wavefunctions are spin-singlets this leaves the 
state invariant. We thus restrict attention to t' < in 
Fig. ^A) . We next turn to various limits in which some 
analytical progress is possible. 



V. SHORT-RANGE RVB LIMIT 

Note that large negative fi with /i < — 4(\t\ + \t'\) is the 
bosonic limitf 1 of the BCS wavefunction, with ip(r — r') 
decaying exponentially in real space. In this limit |0) 
may be viewed as a short-range RVB state. Earlier stud- 
ies of closely related dimer modelsi§*2£i2£ indicate that, 
for short-range singlet bonds there exist four topologi- 
cal sectors on a torus, which may be straightforwardly 
related to sectors with/without visons in the x, y direc- 
tions. In particular, the states with a definite "parity" in 
the dimer model notation correspond to a superposition 
of the states with and without a visoni§*2£i2£. 

Further, two dimer configurations from different sec- 
tors have exponentially small overlap ~ exp(— aL). A 
state is a superposition of many configurations, and as- 
suming that individually small overlaps do not add up 
coherently, the overlap of two states from different sec- 
tors would vanish as L — * 00. 

We have numerically checked (see below) that the over- 
lap vanishes for large L validating the above assumption, 
and implying fractionalization in the short range RVB 
limit. Using the symmetry discussed above, we conclude 
that states in the region with (positive) fi > 4(\t\ + \t'\) 



are also fractionalized, even though ^(r — r') is not obvi- 
ously short-ranged here. 

VI. FURTHER ANALYTICAL INSIGHTS 

We next focus on the curve [i = ^bcs(*') (see 
Fig. ^A)), where ^bcs yields half-filling for the un- 
protected BCS state in the grand canonical representa- 
tion [not fixed particle number representation). This 
unprojected grand canonical state, which we shall de- 
note as \BCS)gc, viewed as the ground state of a 
BCS Bogoliubov-deGennes (BdG) Hamiltonian HsdG: is 
a coherent superposition of number eigenstates sharply 
peaked at the correct mean density. Working with the 
grand canonical wavefunction proves convenient since in 
our discussion below, we will use certain results valid for 
the BdG Hamiltonian in order to infer properties of its 
ground state wavefunction. Note that Gutzwiller pro- 
jecting the grand canonical wavefunction at half-filling 
also picks out the N — L 2 contribution with no double 
occupancy. 

Let us define local SU(2) gauge transformations^ 
U, generated by T+ = c^c) v Tf = Cjl c jh Tf = 

J2a c ]a c ja ~ which mix empty and doubly occupied 
sites as an SU(2) doublet, but act trivially on the 
projected subspace with exactly one particle per site. 
Gutzwiller projection is then equivalent to projection 
onto the SU(2) singlet subspace, and we may write any 
state V\^) — f u U\&), where the integral is over all group 

elements Uq = exp(i J^j Tj ■ 9j). We can thus write 

(U|0)= f {{hc/2e)\U\BCS) GC , (6) 
Ju 

which reduces the problem to computing overlaps of un- 
projected grand canonical states. 

For a nonbipartite lattice (e.g., with t' ^ 0), one cannot 
gauge away the off-diagonal term in the d-wave BCS-BdG 
Hamiltonian by any unitary transformation UHBdGU -1 , 
and the ground state of this transformed Hamiltonian, 
U\BCS)gC: is thus always a SC ground stated (or vortex 
vacuum) for arbitrary U. Then {{hc/2e)\U\BCS) GC = 
for arbitrary U, and thus (V|0) vanishes. Thus nonbi- 
partiteness is a sufficient condition for fractionalization 
in Gutzwiller projected d-wave states^ when /i = /ibcs- 
Note that if one can show this result directly working 
with the wavefunction \BCS) in the fixed number rep- 
resentation (rather than \BCS)gc)i without appeal to 
the BdG Hamiltonian, then this result would be valid for 
general [i and not restricted to // = /ibcs- However, we 
have not found a simple way to show this. 

Let us now turn to the bipartite system (in our case 
t' = 0), where it is well knownS that one can gauge 
away the off-diagonal part of the BdG Hamiltonian 
(provided it does not break time reversal) and trans- 
form the BCS state into a staggered flux state. Now 
({hc/2e)\U\BCS)GC is nonzero for some choice U and one 
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cannot use eq. © to argue that (V\0) vanishes. Note, 
this is not a proof that the overlap remains nonzero in the 
thermodynamic limit since one may worry that this par- 
ticular SU(2) rotation, which converts the d-wave state 
to the staggered flux state, has zero measure in the in- 
tegral and thus still preserve orthogonality. Our numer- 
ical results, presented below, however indicate that the 
overlap does remain nonzero in the thermodynamic limit 
when t' = 0, so that the wavefunction at this special 
point does not describe a Z 2 fractionalized state. 

VII. OVERLAPS WITHIN PURE Z 2 GAUGE 
THEORIES 

Before proceeding to the numerics, we discuss how 
(V|0) is expected to scale with system size in an L x x L y 
system. For ease of presentation, we consider the system 
to be defined on a cylinder as shown in Fig. [21 with peri- 
odic boundary conditions along x. In a unfractionalized 
phase the overlap between the states with and without 
a vison threading the hole of the cylinder will remain 
nonzero in the thermodynamic limit. However, in a frac- 
tionalized phase we argue that on the cylinder the overlap 
should vanish exponentially with L y . This is easy to see 
when the matter fields are gapped (e.g., the short-range 
RVB limit) and may be integrated out to obtain a 2 + 1 
dimensional Z 2 gauge theory. Deep in the deconfining 
phase of the pure gauge theory, a perturbative calcula- 
tion yields (V\0) ~ exp(— L y /£). We obtain this result 
as follows. 

The Hamiltonian for a pure Z 2 gauge theory on the 
square latticed takes the form 

ff* = -*£IK-*£°« ( 7 ) 

where the Pauli matrices a x , a z live on the links of the 
lattice, and □ denotes an elementary plaquette on the 
lattice. This Hamiltonian has an extensive number of 
conserved quantities, namely the set of operators 

G«=II4. (8) 

jei 

where j 6 i denotes the set of sites neighboring site-i, 
all commute with the Hamiltonian, i.e. [G(i), Hz 2 ] = 
for each i. Since G 2 (i) = 1, it is clear that G(i) = ±1 
are the only allowed eigenvalues. For a translationally 
invariant spin liquid insulator with exactly one spin per 
site, it is appropriate to choose G(i) = — 1 at each sit o 5 i 21 
, and this constrained Hamiltonian is referred to as an 
"odd" -gauge theor y 21 ' 35 . Let us represent by P oc id, the 
projection into this subspace, so 

PoM = l[±[l-G(i)]. (9) 

i 

Then, any wavefunction may be projected into the 
"odd" gauge theory sector as P dd\^)- 



(a) (b) 




FIG. 2: Reference configurations to construct the ground 
states for the pure gauge theory with h = (i.e., no vison 
dynamics). (A) The configuration with a z = +1 on all links. 
(B) The configuration with <r z = — 1 on the dark bonds, and 
<j z = +1 on all other links. These two configurations upon 
gauge projection (see text), lead to the two topologically de- 
generate ground states in the gauge theory with h = 0. 



We will now examine the ground states of this Hamilto- 
nian (on a cylinder) deep in the deconfmcd phase, which 
obtains when h/K <C 1. To begin, let us consider the 
extreme limit h = 0. In this case we want to set the flux 
Y[ a o z = 1 on each plaquette to minimize the energy 
of Hz 2 in Eq. (0). One ground state in this limit can 
be achieved by setting a z = 1 on each link as shown in 
Fig. Hl^A), and acting on this reference state with P dd- 
Let us call the resulting state \v = 0). On a cylinder, 
the system has a second, distinct, eigenstate \v — 1); this 
is obtained by gauge projecting the reference state de- 
picted in Fig. |2{B), where a column of horizontal links 
has u z = — 1 such that this reference state also has zero 
flux per plaquette. This state \v = 1) can be alterna- 
tively obtained by acting on the first state \v = 0) with 
the operator 

v+= n (io) 

i£colurnn 

which may be viewed as the operator which creates a Z 2 
vortex (vison) through the hole of the cylinder, since it 
leads to a change in the product of cr z taken around the 
circumference of the cylinder from +1 — > — 1. We may 
loosely refer to it as the one- vison state. 
Clearly, 

(v = 0\v = l) h=0 = 0, (II) 

the overlap between the two states vanishes identically 
even for a finite-size system. Since \V\Hz 2 ] = 0, the 
state \v = 1} is an orthogonal eigenstate degenerate with 
the state \v — 0). This leads to the two-fold topological 
degeneracy on a cylinder. 

Let us now turn on a small nonzero h, so that we are 
still deep in the deconfmed phase h/K <C 1, and calculate 
the perturbative change in each of the two states above. 
Since ha x flips the gauge field on a link from +1 — > — 1, 
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it leads to a superposition within perturbation theory 

\+) + {h/AK)\-) 
^l + h 2 /16K 2 

on each link, and at small h/K the ground state \v = 0)/j 
may be built by taking a reference state which is a di- 
rect product on all links of the above superposition, and 
gauge projecting it using P dd- The second ground state 
\v = l)h is obtained by acting with on \v = Q)h- It 
is straightforward to see that the overlap (v = 0\v = l)h 
vanishes for odd-L y whereas it is nonzero for even L y . 
(This is consistent with the one-vison state carrying a 
non-trivial crystal momentum for odd-L y , which we men- 
tioned earlier.) We can estimate by a direct calculation 
that in fact the overlap for even L y scales as ~ (h/K) L » 
which vanishes exponentially with increasing L y , lead- 
ing to two orthogonal states in the thermodynamic limit. 
Since \V\Hz 2 ] = even for nonzero h, this second 
state is an orthogonal ground state in the thermodynamic 
limit. We again recover the topological degeneracy, and 
see that the no-vison and one-vison states have an overlap 
which vanishes as exp(— L y /^) for sufficiently large sys- 
tem sizes. Similar results have been obtained in Refl3a. 

A simple way to understand the above results is as fol- 
lows. For a two-state system such as a single-particle in 
a double-well potential with a finite barrier, there is a 
nonzero amplitude within perturbation theory for a par- 
ticle localized in one well to be found in the other well 
due to tunneling across the barrier leading to a nonzero 
overlap between states where the particle localized in dif- 
ferent wells. For a collective coordinate such as the vi- 
son, a "string" of length L y , the tunneling between the 
two states, with the vison localized within the cylinder 
or outside the cylinder, proceeds through intermediate 
states where parts of the string lie in the barrier region 
and the overlap is thus exponentially small in L y . 

The behavior of the overlap with gapless matter fields 
is less well studied though we expect it to further suppress 
tunneling and lead to a smaller overlap; numerically (see 
below) we find clear evidence for exponential decay. Note 
that this decay is very different from the decay we might 
expect from the overlap of two general unrelated many- 
body states of a system with fixed density and L x L y sites, 
which would go as ~ exp(—f3L x L y ) for large system sizes. 



VIII. NUMERICAL RESULTS 

We finally discuss numerical results for vison overlaps 
using the variational Monte Carlo method. For computa- 
tional simplicity^ we focus on (V^|V^). We choose t = 1, 
and fix A = 1.25 for all the calculations^ The follow- 
ing results are obtained in regimes where we have given 
analytical arguments above, (i) In the short-range RVB 
regime we find that the wavefunction is Z2 fractional- 
ized; see Fig. E3A) for large negative fi, where the over- 
lap vanishes even on fairly small system sizes, (ii) On the 
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FIG. 3: Size dependence of the overlap {V x \V y )(L) for A = 
1.25 with (A) t' = and (B) t' = -1/4. The different curves 
correspond to overlaps for different values of /j, as indicated. 
An overlap vanishing with increasing system size indicates 
a Z2 fractionalized state, a non- vanishing overlap (e.g., the 
curve in (A) for fi = 0.0) indicates a state which is not Z2 
fractionalized. The lines are fits to the form a(l — tanh((L — 
£*)/£))> which works well for all parameters studied, and from 
which we extract £* . (C) Exponential asymptotic behavior of 
the overlap shown on semi-log plot (see text in Section I VIII 
and Section IVlIll for details). 

curve n = /xbcs we find that nonbipartite (i.e., t' ^ 0) 
wavefunctions are fractionalized as evinced by overlaps 
which vanish with increasing system size (e.g., the curve 
in FigEjB) with t' = -0.25 and // = -0.4). (iii) At the 
special point (t' , /i) = (0,0) with bipartite symmetry, we 
indeed find that the overlaps are nonzero, independent of 
L (see Fig. |3I A)). Thus the vison is not well defined and 
the system is not Z2 fractionalized. Although consistent 
with the analytic arguments, the numerical results are 
important because the arguments were not rigorous for 
(i) and (iii). 

We next evaluate overlaps at t' = 0.0,-0.25,-0.5 for 
— 1.5 < ji < 1.5 which, as we discuss below, covers the 
region of interest for possible spin-liquid insulators in the 
vicinity of high Tc cuprates. We show in Fig.[3]the over- 
lap (V^IVy) as a function of system size L for a range of 
t' and fi. We find that at small L the overlaps are finite 
and then cross over on a scale £* to an asymptotic decay 
exp(— 2L/£). This behavior can be simply described by 
the functional form (V x \V y )(L) = a(l -tanh((L -£*)/£)) 
which appears to fit the data wellS. In particular, this 
form captures the crossover from a constant on small sys- 
tem sizes to an asymptotic exponential form apparent in 
FigHC). 

Deep in the fractionalized regime, we can extract both 
£* and £ and we find £* ~ 2£. However, it is hard to 
access the asymptotic behavior in the region of interest 
around (t' = 0,/i = 0), though we can still reliably ex- 
tract the crossover scale £*. 

We plot the inverse length scale l/£* for t' — 
0, -0.25, -0.5 at different // in Fig. H For t' = 0, we 
see that £* — ► 00 as /i — > 0, fully consistent with the fi- 
nite overlap independent of L on accessible system sizes 



for (t',fi) = (0,0) seen in Fig.^A). It also clear that £* 
remains finite everywhere away from the origin though it 
may become quite large in its vicinity. 

Since fractionalization manifests itself at length scales 
larger than £*, our numerical results strongly suggest 
that all points in the phase diagram of Fig. IHA) are 
fractionalized, except for the origin where £* diverges. 

= (0,0), with its special bipartite symmetry is 
a unfractionalized "singular point" in the space of spin- 
liquid insulators that we study. At this point, the wave- 
function has a power law decay of spin correlations^ 
~ ( — l) :E+i '/|r| Q; with a non-trivial a ~ 1.5, and thus 
appears to be magnetically critical as well. 



i/r 



IX. IMPLICATIONS FOR THE CUPRATES 



Now the question arises — what regime in parameter 
space is relevant for the cuprates? One possibility is that 
there are no fractionalized states in the vicinity of the 
observed phases in real materials. An alternative worth 
exploring in view of the success of RVB wavefunctions 
in understanding the SC stated, is that we take the in- 
sulating limit (hole doping x^O) of V\BCS). Further, 
although V\BCS) is not the ground state of the for the 
half-filled Hubbard model at large-Z7, and does not have 
long-range antiferromagnetic order, it is known to be an 
energetically competitive candidate statei^S^ 

Optimizing variational parameters /i, A for the large- 
ly Hubbard model (U= 12, t'= -1/4), we find A ~ 1.25 
while —0.3 <■ fi < 0.3 as shown in Fig. H^B). For this 
region (the shaded ellipse in Fig.^A)) we conclude from 
Fig. ElB) that £* > 25 lattice spacings. 

We now convert this to an estimate of the energy scale 
E v below which fractionalization is apparent. We ex- 
pect E v — aJ/(£_*) z < aJ/£* which vanishes at the 
"singular point" with bipartite symmetry. Here J is the 
nearest-neighbor superexchange, a = 0(1) is a dimen- 
sionless constant and the dynamical exponent z > 1. For 
a = 1 and J = 1200if, the estimated E v < 50K. More 
concretely, proximity to the unfractionalized 'bipartite' 
point (t' = 0, (j, = 0) can lead to a very small vison gap 
for fractionalized RVB states. Recently, based on an inge- 
nieous proposal of Senthil and Fisher 3 ™, a flux trapping 
experiment was been carried out to detect the vison 
in highly underdoped cuprates. However, these experi- 
ments did not see the vison, which led to an upper bound 
on the vison gap < 150-Rf . Our estimates are consistent 
with this bound, and suggest that Z2 fractionalization, 
even if present in the cuprates, would likely be apparent 
only at very low T and cannot play a role in pseudogap 
anomalies. 
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FIG. 4: Dependence of l/£* on fi for (A) t' = 0.0, (B) t' = 
—0.25, (C) t' = —0.5. In panel (A), l/£* extrapolates to zero 
on approaching (t' , [i) = (0,0) indicating a unfractionalized 
"singular point" at the origin. For t' — —0.25, —0.5 (B,C) the 
length scale £* appears finite everywhere. 

X. APPLICATION TO NUMERICAL STUDIES 
OF HUBBARD-LIKE MODELS 



At this stage, we would like to emphasize the most im- 
portant general outcome of our numerical results which 
goes beyond the specific conclusions regarding the par- 
ticular wavefunction we have studied, and how it can be 
applied to numerical studies of model system; 
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test for Z2 fractionalization. If a system is in a Z2 frac- 
tionalized insulating phase, we have shown that the over- 
lap of the wavefunctions with periodic and antipcriodic 
boundary conditions imposed on the electrons must van- 
ish exponentially as exp(— L/£) on an L x L system. The 
length scale beyond which fractionalization is apparent, 
£*, can be deduced from studying the finite size scaling 
of this overlap. These two results are useful in detect- 
ing topological order when the full low lying spectrum of 
a microscopic Hamiltonian is not accessible so that the 
topological degeneracies are not obvious even though one 
can work with large system sizes. 

Consider, for example, Lanczos studies which improve 
upon a definite trial wavefunction for any Hamiltonian 
using a few Lanczos iterations^, or other Monte Carlo 
studies^ which also use wavefunctions. Let us imagine 
working on a cylinder (or torus), and starting with two 
different wavefunctions which differ in the boundary con- 
dition on the electrons. In a Z2 spin liquid insulator, this 
should provide us with the two topologically degenerate 
ground states. However, if the Hamiltonian has a conven- 
tional ground state however, we expect the two resulting 
states to be identical. 

Note that, if one has small system sizes, the two initial 
states may have a large overlap even in a Z2 spin liquid, 
and then would lead to a single ground state. However, 
on a large enough system, we expect the initial overlap 
would be exponentially small in a Z2 spin liquid, so that 
one may recover the different topological sectors in this 
manner. Another check, apart from a vanishing overlap 
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of the different final states, is that the energy difference 
between them should vanish with increasing system size. 
If the spinon excitations are gapped the energy difference 
is expected to decay exponentially in L, but it would 
have a power law decay in the presence of gapless spinon 
excitations. 



should be applied to these systems. If they are shown 
to have topological order consistent with that expected 
for a Zi spin liquid, they may provide us with the first 
examples of Zi fractionalization in microscopic models 
with full spin rotational symmetry. 



XI. CONCLUDING REMARKS 

In this paper, we have analyzed the Gutzwiller pro- 
jected d-wave BCS state at half-filling, and shown from 
analytic arguments and numerical results that it gener- 
ically describes a Z^ fractionalized spin liquid except at 
a special 'bipartite' point where the vison is no longer 
well defined. We have pointed out the significance of this 
result for the cuprate superconductors. As mentioned 
in the introduction, there have been many recent exam- 
ples of spin-liquid states reported in numerical studies 
of SU{2) spin models^ and Hubbard-like modelo 14 ! 15 . 
These models have been identified as having spin liquid 
ground states based on the fact that the ground states do 
not appear to possess any simple broken symmetry pat- 
terns. The ideas and methods developed in this paper 
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